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Quantum heat engines (QHE) are thermal machines 
where the working substance is quantum. In the extreme 
case the working medium can be a single particle or a few 
level quantum system. The study of QHE has shown a re¬ 
markable similarity with the standard thermodynamical 
models, thus raising the issue what is quantum in quan¬ 
tum thermodynamics. Our main result is thermodynam¬ 
ical equivalence of all engine type in the quantum regime 
of small action. They have the same power, the same 
heat, the same efficiency, and they even have the same 
relaxation rates and relaxation modes. Eurthermore, it 
is shown that QHE have quantum-thermodynamic signa¬ 
ture, i.e thermodynamic measurements can confirm the 
presence of quantum coherence in the device. The coher¬ 
ent work extraction mechanism enables power outputs 
that greatly exceed the power of stochastic (dephased) 
engines. 


I. INTRODUCTION 

Thermodynamics emerged as a practical theory for 
evaluating the performance of steam engines. Since then 
the theory proliferated and utilized in countless systems 
and applications. Eventually, thermodynamics became 
one of the pillars of theoretical physics. Amazingly, it 
survived great scientific revolutions such as quantum me¬ 
chanics and general relativity. In fact, thermodynamics 
played a crucial role in the development of these theories. 
Black body radiation led Planck and Einstein to intro¬ 
duce energy quantization, and the second law of thermo¬ 
dynamics led Bekenstein to discover the relation between 
the event horizon area and black hole entropy and tem¬ 
perature. 

The Carnot efficiency is a manifestation of the sec¬ 
ond law in heat engines, and it is universally valid. In 
addition all reversible engines must operate at Carnot 
efficiency. Though very profound these principles have 
limited practical value. Eirst, real engines produce fi¬ 
nite non zero power and therefore cannot be reversible. 
Second, the performance of real engine is more severely 
limited by heat leaks, friction and heat transport. This 
led to the study of efficiency at maximal power mm and 
finite time thermodynamics mill- In finite time thermo¬ 
dynamics time was introduced classically using empirical 
heat transport models. Eor microscopic systems the nat¬ 
ural avenue to introduce time dynamics is via quantum 
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mechanics of open systems. The full quantum dynamics 
may lead to new mechanism of extracting work or cool¬ 
ing. Alternatively, it may lead to new microscopic heat 
leaks and friction-like mechanisms. 

Quantum thermodynamics is the study of thermody¬ 
namic quantities such as temperature, heat, work, and 
entropy in microscopic quantum systems or even for a sin¬ 
gle particle. This includes dynamical analysis of engines 
and refrigerators in the quantum regime EMIT], theoreti¬ 
cal frameworks that takes into account single shot events 
[2811^ . and the study of thermalization mechanisms [30^ 
|32]. 

It was natural to expect that in the quantum regime 
new thermodynamic effects will surface. However, quan¬ 
tum thermodynamic systems (even with a single particle) 
show a remarkable similarity to macroscopic system de¬ 
scribed by classical thermodynamic. When the baths are 
thermal the Carnot efficiency limit is equally applicable 
for a small quantum system [5j [33] . Even classical ffac¬ 
tuation theorems hold without any alteration [511155] . 

Is there really nothing new and profound in the ther¬ 
modynamics of small quantum system? In this work, we 
present thermodynamic behavior that is purely quantum 
in its essence and has no classical counterpart. 

Recently some progress on the role of quantum coher¬ 
ence in quantum thermodynamics has been made [371158]. 
In addition, quantum coherence has been shown to quan¬ 
titatively affect the performance of heat machines 
m]. In this work we associate coherence with a specific 
thermodynamic effect and relate it to a thermodynamic 
work extraction mechanism. 

Heat engines can be classified by their different 
scheduling of the interactions with the baths and the 
work repository. These types include the four-stroke, 
two-stroke and the continuous engines (these engine 
types will be described in more detail later on). The 
choice of engine type is usually guided by convenience of 
analysis or ease of implementation. Nevertheless, from 
a theoretical point of view, the fundamental differences 
or similarities between the various engine types are still 
uncharted. This is particularly true in the microscopic 
quantum regime. Eor brevity we discuss engines but all 
our results are equally applicable to other heat machines 
such as refrigerators and heaters. 

Our first result is that in the limit of small engine 
action (weak thermalization, and a weak driving field), 
all three engine types are thermodynamically equivalent. 
The equivalence holds also for transients and for states 
that are very far from thermal equilibrium. On top of 
providing a thermodynamic unification limit for the var¬ 
ious engine types, this finding also establishes a connec- 
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tion to quantum mechanics as it crucially depends on 
phase coherence and quantum interference. In particu¬ 
lar, the validity regime of the equivalence is expressed in 
terms of h. 

Our second result concerns quant um-t her mo dynamic 
signatures. Let us define a quantum signature as a signal 
extracted from measurements that unambiguously indi¬ 
cates the presence of quantum effects (e.g. entanglement 
or interference). The Bell inequality for the EPR ex¬ 
periment is a good example. A quantum-thermodynamic 
signature is a quantum signature obtained from measur¬ 
ing thermodynamic quantities. We show that it is pos¬ 
sible to set an upper bound on the work output of a 
stochastic, coherence-free engine. Any engine that sur¬ 
passes this bound must have some level of coherence. 
Hence, work exceeding the stochastic bound constitutes 
a quantum-thermodynamic signature. Furthermore, we 
distinguish between a coherent work extraction mecha¬ 
nism and a stochastic work extraction mechanism. This 
explains why in the equivalence regime, coherent engines 
produce significantly more power compared to the corre¬ 
sponding stochastic engine. 

The equivalence derivation is based on three ingredi¬ 
ents. First, we introduce a multilevel embedding frame¬ 
work that enables the analysis of all three types of engines 
in the same physical setup. Next, a “norm action” small¬ 
ness parameter, s, is defined for engines using Liouville 
space. The third ingredient is the symmetric rearrange¬ 
ment theorem that is used to show why all three engine 
types have the same thermodynamic properties despite 
the fact that they exhibit very different density matrix 
dynamics. 

In section two we describe the main engine types, and 
introduce the multilevel embedding framework. Next, 
in section three the multilevel embedding and the sym¬ 
metric rearrangement theorem are used to show the var¬ 
ious equivalence relation of different engine types. After 
discussing the two fundamental work extraction mech¬ 
anisms, in section four we present and study the over- 
thermalization effect in coherent heat engines. In section 
five we show a quantum thermodynamic signature that 
separates quantum engines from stochastic engines. Fi¬ 
nally, in section six we conclude and discuss extensions 
and future prospects. 


II. HEAT ENGINES TYPES AND THE 
MULTILEVEL EMBEDDING SCHEME 

Heat engines are either discrete such as the two stroke 
engines or the four stroke Otto engine. Heat engines can 
also operate continuously such as in turbines. Quantum 
analogues of all these engine types have been studied [73]. 
Here we present a theoretical framework where all three 
type of engines can be embedded in a unified physical 
framework. This framework termed “multilevel embed¬ 
ding” is an essential ingredient in our theory as it enables 
a meaningful comparison between different engine types. 


A. Heat and work 

A heat engine is a device that uses at least two thermal 
baths in different temperatures to extract work. Loosely 
speaking, work is the transfer of energy to a single de¬ 
gree of freedom. For example, increasing the excitation 
number of an oscillator, increasing the photon number 
in a specific optical mode (lasing) or increasing the ki¬ 
netic energy in a single predefined direction. “Battery” 
or “flywheel” are terms often used in this context of work 
storage mi m- We shall use the more general term 
“work repository”. Heat, on the other hand, is energy 
spread over multiple degrees of freedom. Close to equi¬ 
librium and in quasistatic process, the heat is related 
to the temperature and to the effective number of de¬ 
grees of freedom (entropy S) via the well known relation 
dQ = TdS. 

In the elementary quantum heat engines the working 
substance comprises of a single particle (or few at the 
most). Thus the working substance cannot reach equilib¬ 
rium on its own. Furthermore, excluding few non-generic 
cases it is not possible to assign an equation of state that 
establishes a relation between thermodynamic quantities 
when the substance is in equilibrium. Nevertheless, QHE 
satisfy the second law and therefore are also bounded by 
Carnot efficiency limit [5]. 

Work strokes are characterized by zero contact with 
the baths and an inherently time dependent Hamiltonian. 
The unitary evolution generated by this Hamiltonian can 
change the energy of the system. On the other hand, the 
von Neumann entropy and the purity, remain fixed (uni¬ 
tary evolution at this stage). Hence the energy change of 
the system in this case constitutes pure work. The sys¬ 
tem’s energy change is actually an energy exchange with 
the work repository. 

When the system is coupled to a thermal bath and the 
Hamiltonian is fixed in time, the bath can change the 
populations of the energy levels. In steady state the sys¬ 
tem reaches a Gibbs state where the density matrix has 
no coherences in the energy basis and the population of 

the levels is given by: ^ / 2n=i ^ where 

N is the number of levels and 'h' stands for 'c' (cold) or 
'h' (hot). In physical models where the system thermal- 
izes via collision with bath particles a full thermalization 
can be achieved in finite time [44II48] . However, it is not 
necessary that the baths will bring the system close to 
a Gibbs state for the proper operation of the engine. In 
particular, maximal efficiency (e.g. in Otto engines) can 
be achieved without full thermalization. Maximal power 
(work per cycle time) is also associated with partial ther¬ 
malization [TJ|49]. The definitive property of a thermal 
bath is its aspiration to bring the system to a predefined 
temperature regardless of the initial state of the system. 
The evolution in this stage does not conserve the eigen¬ 
values of the density matrix of the system, and therefore 
not only energy but entropy as well is exchanged with 
the bath. Therefore, the energy exchange in this stage is 
considered as heat. 
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In contrast to definitions of heat and work that are 
based on the derivative of the internal energy [5l[50l[5l], 
our definitions are obtained by energy balance when cou¬ 
pling only one element (bath or external field) at a time. 
As we shall see later on, in some engine types several 
agents change the internal energy simultaneously. Even 
in this case, this point of view of heat and work will still 
be useful for obtaining consistent and physical definitions 
of heat and work. 


B. The three engine types 

There are three core engine types that operate with two 
thermal baths: four-stroke engine, two-stroke engine, and 
a continuous engine. A stroke is a time segment where 
a certain operation takes place, for example thermaliza- 
tion or work extraction. Each stroke is a CP map and 
therefore the one-cycle evolution operator of the engine 
is also a CP map (since it is a product of CP maps). Eor 
the extraction of work it is imperative that some of the 
stroke propagators do not commute [5^ . 

Otto engines and Carnot engines are examples of four 
stroke engines. The simplest quantum four-stroke en¬ 
gine is the two-level Otto engine shown in Eig. la. In 
the first stroke only the cold bath is connected. Thus, 
the internal energy changes are associated with heat ex¬ 
change with the cold bath. The expansion and compres¬ 
sion of the levels are fully described by a time-dependent 
Hamiltonian of the form H(t) = f{t)az (the baths are 
disconnected at this stage). In the second stroke, work 
is consumed in order to expand the levels, and in the 
fourth stroke work is produced when levels revert to their 
original values. There is a net work extraction since the 
populations in stages II and IV are different. As we shall 
see later on, other unitary operations are more relevant 
for quantum equivalence of heat engines. Nevertheless, 
this particular operation resembles the classical expan¬ 
sion and compression of classical engines. The work is 
the energy exchanged with the system during the unitary 
stages: W = Wn + Wiv = ((^ 3 )-(^ 2 ))+ ((^ 5 )-(^ 4 )). 
We will consider only energy expectation values for two 
main reasons. Eirst, investigations of work fluctuations 
revealed that quantum heat engines follow classical fluc¬ 
tuation laws [ 53 ] and we search for quantum signatures in 
heat engines. The second reason is that in our view the 
engine should not be measured during operation. The 
measurement protocol used in quantum fluctuation the¬ 
orems [smaisi], eliminates the density matrix coher¬ 
ences. These coherences have a critical component in 
the equivalence and quantum signature we study in this 
paper. Thus, although we frequently calculate work per 
cycle, the measured quantity is the cumulative work and 
it is measured only at the end of the process. The av¬ 
eraged quantities are obtained by repeating the full ex¬ 
periment many times. Engines are designed to perform a 
task and we assume that this completed task is the sub¬ 
ject of measurement. The engine internal state are not 



Figure 1: (a) A two-level scheme of a four-stroke engine, (b) 
Two-particle scheme of a two stroke engine, (c) A three-level 
scheme of a continuous engine. 


measured. 

The heat per cycle taken from the cold bath is Qc = 
{E 2 ) — {El) and the heat taken from the hot bath is 
Qh = (E 4 ) — (E 3 ). In steady state the average energy of 
the system returns to its initial value after one cycleso 
that (E 5 ) = (^i)- Erom this it follows immediately that 
Qc + Qh + VE = 0 , i.e. the first law of thermodynamics is 
obeyed. There is no instantaneous energy conservation 
of internal energy, as energy may be temporarily stored 
in the interaction field or in the work repository. 

In the two-stroke engine shown in Eig lb the engine 
consists of two parts (e.g. two qubits) [5l]. One part 
may couple only to the hot bath and the other may cou¬ 
ple only to the cold bath. In the first stroke both parts 
interact with their bath (but do not necessarily reach 
equilibrium). In the second unitary stroke the two engine 
parts are disconnected from the baths and are coupled 
to each other. They undergo a mutual unitary evolution 
and work is extracted in the process. 

In the continuous engine shown in Eig. Ic the two 
baths and the external interaction field are connected 
continuously. Eor example, in the three level laser sys¬ 
tem shown in Eig Ic the laser light represented by Pwit) 
generates stimulated emission that extracts work from 
the system. This system was first studied in thermo¬ 
dynamics context in m, while a more comprehensive 
analysis of the system was given in m- It is imperative 
that the external field is time-dependent. If it is time- 
independent the problem becomes a pure heat transport 
problem where Qh = —Qc 9 ^ 0. In heat transport the 
interaction field merely “dresses” the level so that the 
baths see a slightly modified system. The Lindblad gen¬ 
erators are modified accordingly and heat flows without 
extracting or consuming work m- Variations on these 
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engine types may emerge due to realization constraints. 
For example, in the two stroke engine the baths may be 
continuously connected. This variation and others can 
still be analyzed using the tools presented in this paper. 

C. Efficiency vs. work and heat 

Since the early days of Carnot, efficiency received con¬ 
siderable attention for two main reasons. First, this 
quantity is of great interest from both theoretical and 
practical points of view. The second reason is that effi¬ 
ciency satisfies a universal bound that is independent of 
the engine details. The Carnot efficiency bound is a man¬ 
ifestation of the second law of thermodynamics. Indeed, 
for Markovian bath dynamics it was shown that quantum 
heat engines cannot exceed the Carnot efficiency [5]. Re¬ 
cently, a more general approach based on a fluctuation 
theorem for QHE showed that the Carnot bound still 
holds for quantum engines m Studies in which higher 
than Carnot efficiency are reported gD], are interesting 
but they use non-thermal baths and therefore not surpris¬ 
ingly deviate from results derived in the thermodynamic 
framework that deals with thermal baths. For example, 
an electric engine is not limited to Carnot efficiency since 
its power source is not thermal. Although, the present 
work has an impact on efficiency as well, we focus on work 
and heat separately in order to unravel quantum effects. 
As will be exemplified later, in some elementary cases 
these quantum effects do not influence the efficiency. 

D. Bath description and Liouville space 

The dynamics of the working fluid (system) interact¬ 
ing with the heat baths is described by Lindblad-Gorini- 
Kossakowski-Sudarshan (LGKS) master equation for the 
density matrix [581I6Q] : 

dtp = L{p) = -i[H„p]+J2^kpAl-^AlAkp-^pAlAk, 

k 

( 1 ) 

where the operators depend on the temperature, re¬ 
laxation time of the bath, system bath coupling, and 
also on the system Hamiltonian Hg [58]. This form al¬ 
ready encapsulates within the Markovian assumption of 
no memory. The justification for these equations arises 
from a “microscopic derivation” in the weak system bath 
coupling limit [61]. In this derivation a weak interaction 
field couples the system of interest to a large system (the 
bath) with temperature T. This interaction brings the 
system into a Gibbs state at temperature T. The Lind- 
blad thermalization operators used for the baths are 
described in the next section 

Equation 0 is a linear equation so it can always be 
rearranged into a vector equation. Given an index map¬ 
ping PnxN \p)ixN^ Lindblad equation now reads: 

idt \p) = {ns + C) Ip) , (2) 


I 



Figure 2: In the standard two-level Otto engine there are two- 
level Eg^e (purple) that change in time to Eg In the mul¬ 
tilevel embedding framework, the levels (E 1 - 4 ) are fixed in 
time (black dashed lines) but a time dependent field (tt pulse, 
swap operation) transfer the population (green) to the other 
levels. For a swap operation the two schemes lead to the same 
final state and therefore are associate with the same work. 
Nonetheless, the multilevel scheme is more general since for 
weaker unitary transformation (instead of the tt pulse), co¬ 
herences are generated. We show that this type of coherences 
can significantly boost the power output of the engine. 

where ELs is an Hermitian x matrix that originates 
from Hs, and £ is a non Hermitian x matrix 
that originates from the Lindblad evolution generators 
Ak. This extended space is called Liouville space [ 6 ^ . 
In this paper we will use calligraphic letters to describe 
operators in Liouville space and ordinary letters for op¬ 
erators in Hilbert space. For states, however, \A) will 
denote a vector in Liouville space formed from Anxn 
by “vec-ing” A into a column in the same procedure p is 
converted into \p). A short review of Liouville space and 
some of its properties is given in appendix H. 

While the Lindblad description works very well for suf¬ 
ficiently long times it fails for very short times where 
some of the approximation breaks down. In scales where 
the bath still has a memory of the system’s past states, 
the semi group property of the Lindblad equation no 
longer holds: \p{tEt')) 7 ^ ^ \p{t')). This 

will set a cutoff limit for the validity of the engine types 
equivalence in the Markovian approximation. 

Next we introduce the multilevel embedding scheme 
that enables us to discuss various heat engines in the 
same physical setup. 


E. Multilevel embedding 

Let the working substance of the quantum engine be an 
A/'-level system. These levels are fixed in time (i.e. they 
do not change as in Fig. la). For simplicity, the levels 
are non degenerate. We divide the energy levels into a 
cold manifold and a hot manifold. During the operation 
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Figure 3: Representation of the three types of engines in the 
multilevel embedding framework. In this scheme the different 
engine types differ only in the order of coupling to the baths 
and work repository. Since the interactions and energy levels 
are the same for all engine types, a meaningful comparison of 
performance becomes possible. 


of the engine the levels in the cold manifold interact only 
with the cold bath, and the hot manifold interact only 
with the hot bath. Each thermal coupling can be turned 
on and off as a function of time but the aliasing of a level 
to a manifold does not change in time. 

If the manifolds do not overlap the hot and cold ther¬ 
mal operations commute and they can be applied at the 
same time or one after the other. The end result will 
be the same. Nevertheless, our scheme also includes the 
possibility that one level appears in both manifold. This 
is the case for the three-level continuous engine shown 
in Fig Ic. For simplicity, we exclude the possibility of 
more than one mutual level. If there are two or more 
overlapping levels there is an inevitable heat transport 
in steady state from the hot bath to the cold bath even 
in the absence of an external field that extract work. In 
the context of heat engines this can be interpreted as heat 
leak. This “no field - no transport” condition holds for 
many engines studied in the literature. Nonetheless, this 
condition is not a necessary condition for the validity of 
our results. 

This manifold division seems sensible for the continu¬ 
ous engine and even for the two stroke engine in Fig lb, 
but how can it be applied to the four stroke engine shown 
in Fig. la? The two levels interact with both baths and 
also change their energy value in time contrary to the as¬ 
sumption of fixed energy levels. Nevertheless, this engine 
is also incorporated in the multilevel embedding frame¬ 
work. Instead of two-levels as in Fig. la consider the 
four-level system shown in the dashed green lines in Fig 
2 . 

Initially, only levels 2 and 3 are populated and coupled 
to the cold bath (2 & 3 are in the cold manifold). In the 
unitary stage an interaction Hamiltonian Hgwap gener¬ 
ates a full swap of populations and coherence according 
to the rule 1 2, 3 4. Now levels 1 and 4 are pop¬ 

ulated and 2 and 3 are empty. Therefore, this system 
fully simulates the expanding level engine shown in Fig. 
la. At the same time this system satisfies the separation 


into well defined time-independent manifolds as defined 
in the multilevel embedding scheme. 

The full swap used to embed the traditional four-stroke 
Otto engine, is not mandatory and other unitary oper¬ 
ations can be applied. This extension of the four-stroke 
scheme is critical for our work since the equivalence of 
engines appear when the unitary operation is fairly close 
to the identity transformation. 

Figure 3 shows how the three engine types are rep¬ 
resented in the multilevel embedding scheme. The ad¬ 
vantage of the multilevel scheme now becomes clear. All 
three engine types can be described in the same physical 
system with the same baths and the same coupling to 
external fields (work extraction). The engine types differ 
only in the order of the coupling to the baths and to the 
work repository. While the thermal operations commute 
if the manifolds do not overlap, the unitary operation 
never commutes with the thermal strokes. 

On the right of Fig. 3, we plotted a “brick” diagram for 
the evolution operator. Black stands for unitary trans¬ 
formation generated by some external field, while blue 
and red stand for hot and cold thermal coupling, respec¬ 
tively. When the bricks are on top of each other it means 
that they operate simultaneously. Now we are in position 
to derive the first main results of this paper: the ther¬ 
modynamic equivalence of engine types in the quantum 
regime. 

III. CONTINUOUS AND STROKE ENGINE 
EQUIVALENCE 

We first discuss the equivalence of continuous and four- 
stroke engines. Nevertheless, all the argument are valid 
for the two-stroke engines as well, as explained later on. 
Although our results are not limited to a specific en¬ 
gine model it will be useful to consider the simple engine 
shown in Fig. 4 . We will use this model to highlight 
a few points and also for numerical simulations. The 
Hamiltonian part of the system is: 

Ho ^ cos{ujt)H^, (3) 

where Ho = |1) (1| - ^ |2) (2| + ^ |3) (3| + 

|4) (4|, Hu, = €{t) |1) (2| + e{t) |3) (4| + h.c. and to = 

AEh — AEc 

2 

The driving frequency, that couples the system to the 
work repository is in resonance with the top and bottom 
energy gaps. The specific partitioning into hot and cold 
manifolds was chosen so that only one frequency (e.g. a 
single laser) is needed for implementing the system in¬ 
stead of two. 

We assume that the Rabi frequency of the drive e is 
smaller than the decay time scale of the baths e <C 7 c 5 7/i- 
Under this assumption the dressing effect of the driving 
field on the system-bath interaction can be ignored. It 
is justified, then, to use “local” Lindblad operators ob¬ 
tained in the absence of a driving field [saisi]. For plot¬ 
ting purposes (reasonable duty cycle) in the numerical 
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Figure 4: Illustration of the engine used in the numerical 
simulation. By changing the time order of the coupling to 
and to thermal baths, all three type of engines can realized in 
the model. 

examples we will often use e = 7c = 7/^. While this poses 
no problem for stroke engines realizations, for experimen¬ 
tal demonstration of equivalence with continuous engines 
one has to increase the duty cycle so that e ^ 7c:7/i- 
That is, the unitary stage should be made longer but 
with weaker driving field. 

The Lindblad equation is given by 0 with the Hamil¬ 
tonian ^ and with the following Lindblad operators: 

Ai = v^e |4)(1|, 

^2 = x/^|l)(4|, 

AEc: 

^3 = 2Te |3)(2|, 

^4 = y^|2)(3|. 


Note that Ch^c were not modified by the transformation 
to the rotating system since [Ch^ci"l~Lo] = 0 in the mi¬ 
croscopic derivation [75]. Now that we have established a 
regime of validity and the super Hamiltonian that governs 
the system, we can turn to the task of transforming from 
one engine type to other types and study what properties 
change in this transformation. The engine type transfor¬ 
mation is based on the Strang decomposition m for two 
non commuting operators A and B (the operators may 
not be Hermitian): 

^{A+B)dt _ ^\Adt^Bdt^\Adt _j_ ^\Adt^Bdt^\Adt 

(5) 

where s = (||7l|| + ||S||)(it must be small for the expan¬ 
sion to be valid. ||7l|| is the spectral norm (or opera¬ 
tor norm) of 7l, and it is the largest singular value of 
71, ||7l|| = max\/ eigfTlTltf [65]. For Hermitian operators 
with eigenvalues the spectral norm is max{\\j^^i\). 
In Appendix I we derive the condition s | for the va¬ 
lidity of 0- We will use the symbol = to denote equality 
with correction 0{s^). Let the evolution operator of the 
continuous engine over the chosen cycle time Tcyc = ^'mrd 
be: 

^cont _ (0) 

By first splitting Cc and then splitting Ch we get: 

^4 stroke _ g-i(3£e) ^ ^ g-i(3£^) ^ 


In all the numerical simulation we use /S.Eh = 4, AF^c = 
1, T/i = 5, Tc = 1. The interaction with the baths or 
with work repository can be turned on and off at will. 

Starting with the continuous engine, we choose a unit 
cell that contains exactly 6m (m is an integer) complete 
cycles of the drive {rd = 2 t: juj) so that Tcyc = The 

difference between the engine cycle time and the cycles 
of the external drive, will become clear in stroke engines 
(also the factor of six will be clarified). 

For the validity of the secular approximation used in 
the Lindblad microscopic derivation [58], the evolution 
time scale must satisfy r ^ min{AEh ae ) • Therefore 
m must satisfy m ^ min{AEh ae ) • hlote that if the 
Lindblad description is obtained from a different phys¬ 
ical mechanism (e.g. thermalizing collisions) then this 
condition is not required. 

Next we transform to the interaction picture (denoted 
by tilde) using the transformation U = and per¬ 

form the rotating wave approximation (RWA) by drop¬ 
ping terms oscillating in frequency of 2uj. For the RWA 
to be valid the amplitude of the field must satisfy e uj. 
The resulting Liouville space super Hamiltonian is: 

n = Cc + Ch + \n^. (4) 


Note that the system is periodic so the first and last 
stages are two parts of the same thermal stroke. Con¬ 
sequently 0 describes an evolution operator of a four 
stroke engine, where the unit cell is symmetric. This 
splitting is illustrated in Fig. 5a and Fig. 5b. There are 
two thermal strokes and two work strokes that together 
constitute an evolution operator that describes a four- 
stroke engine. The cumulative evolution time as written 
above is (m + m + 2m + m + m)rc^ = 6mrd = Tcyc- Yet, to 
maintain the same cycle time as chosen for the continuous 
engine, the coupling to the baths and field were multi¬ 
plied by three. In this four-stroke engine each thermal 
or work stroke operates in total only a third of the cy¬ 
cle time compared to the continuous engine. Hence, the 
coupling must be three times larger in order to generate 
the same evolution. 


By virtue of the Strang decomposition IC^ stroke ^ 
^cont ^ ^ 2 . The action parameter s of the en¬ 


gine is defined as 5 = 


cyc/2 
yc/^ 


‘ cyc' 


u dt = ilwu^w + 

Although we used h = 1 units, note 
that s has dimensions of h and therefore the relation 
^4 stroke ^ ^cont ^heu the engine action is 

small compared to h. This first appearance of a quan¬ 
tum scale will be discussed later on. 
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A. Dynamical aspect of the equivalence 

Before discussing the thermodynamics properties of 
the engine we point out that stroke ^ ^cont 
immediate important consequences. First both engines 
have the same steady state solution over one cycle \ps): 

iC^ “(r,,,) \p,) - \Ps) = \Ps), (8) 

{Cc ^ Ch ^-^1-Lw)\ps) = 0 - ( 9 ) 

At time instances that are not an integer multiples Tcyc^ 
the states of the engines will differ significantly (0(5^)) 
since ^4 stroke ^ < T^yc)- That is, the 

engines are still significantly different from each other. 
The second consequence is that the two engines have the 
same transient modes as well. When monitored at mul¬ 
tiple of Tcyc both engines will have the same relaxation 
dynamics to steady state if they started from the same 
initial condition. In the reminder of the paper, when the 
evolution operator is written without a time tag, it means 
we consider the evolution operator of a complete cycle. 


B. Thermodynamic aspect of the equivalence 


The equivalence of the one cycle evolution operators 
of the two engines does not immediately imply that the 
engines are thermodynamically equivalent. Generally, in 
stroke engines the heat and work depend on the dynam¬ 
ics of the state inside the cycle which is very different 
( 0 ( 5 ^)) from the constant state of the continuous engine. 
However, in this section we show that all thermodynam¬ 
ics properties are equivalent in both engines up to 0{s^) 
corrections, similarly to the evolution operator. We start 
by evaluating the work and heat in the continuous en¬ 
gine. By considering infinitesimal time elements where 
>Cc, Ch and Hw operate separately, one obtains that the 
heat and work currents are[76] jc{h) = {Hq |>^c(/i)| Ps{t)) 
and = (i7o | \ Ps{t))' la Itie continuous engine the 

steady state satisfies \ps{t)) = \ps) so the total heat and 
work in steady state in one cycle are: 


^cont 

Q cont 

clh) 




Ps 


'^cyc-) 


(^Hq |>^c(/i) I Ps) 


'T'cyc' 


( 10 ) 

( 11 ) 


Figure 5: (color online) Illustration of the splitting of the evo¬ 
lution operator of the continuous engine (a), into a four-stroke 
engine (b), a two-stroke engine (c), and a two-field four-stroke 
engine (d). The horizontal axis corresponds to time as before. 
The size of the brick corresponds to the strength of coupling 
to the work repository or to the baths. The symmetric rear¬ 
rangement theorem ensures that in the limit of small action, 
any rearrangement that is symmetric with respect to the cen¬ 
ter, and conserves the area of each color does not change the 
total power and heat. 


These quantities should be compared to the work and 
heat in the four-stroke engine. Instead of carrying 
out the explicit calculation for this specific four-stroke 
splitting we use the symmetric rearrangement theorem 
(SRT) derived in Appendix III. Symmetric rearrange¬ 
ment of a Hamiltonian means is a change in the or¬ 
der of couplings e(t), 7 c(t), 7 /i(t) that satisfies f e{t)dt = 
const, f = const, f = const and e(t) = 

e(-i),7c(i) = 

']-[!¥ stroke aud auy other super Hamiltonian obtained 
using the Strang splitting of the continuous engine, are 













Figure 6: The equivalence of heat engine types in transient 
evolution when the engine action is small compared to h. (a) 
The cumulative power transferred to the work repository is 
plotted as a function of time. All engines start in the excited 
state |4), which is very far from the steady state of the system. 
At complete engine cycles (vertical lines) the power in all 
engines is the same, (b) Once the action is increased (here 
the field e was increased), the equivalence no longer holds. 


examples of symmetric rearrangements. The SRT ex¬ 
ploits the symmetry of the Hamiltonian to show that 
symmetric rearrangement change heat and work only in 
0{s^). In the Appendix III we show that 


1^4 stroke j^cont 

Q 4 stroke /<^cont 

c(h) — ^clh) • 


( 12 ) 

(13) 


Thus, we conclude that up to corrections, the engines 
are thermodynamically equivalent. When 5 ^ 1 . work, 
power heat and efficiency converge to the same value for 
all engine types . Clearly, inside the cycle the work and 
heat in the two engine are significantly different ( 0 ( 5 ^)) 
but after a complete cycle they become equivalent. The 
symmetry makes this equivalent more accurate as it holds 
up to (rather s^). Interestingly the work done in the 
first half of the cycle is + 0(<s^). However, when 

the contribution of the other half is added added the 
correction cancels out and ( 12 ) is obtained (see Appendix 
HI). 

We emphasize that the SRT and its implications 
111 ),® are valid for transients and for any initial state 
- not just for steady state operation. In Fig. 6 a we 
show the cumulative work as a function of time for a 
four-stroke engine and a continuous engine. The vertical 
lines indicate complete cycle of the four-stroke engine. 


In addition to the parameter common to all examples 
specified before, we used e = = "fh = 10 “^ and the 

equivalence of work at the vertical lines is apparent. In 
Fig. 6 b the field and thermal coupling were increased 
to e = 7c = 7/i = 5 X 10“^ . Now the engines perform 
differently even at the end of each cycle. This example is 
a somewhat extreme situation where the system changes 
quite rapidly (consequence of the initial state we chose). 
In other cases, such as steady state operation, the equiv¬ 
alence can be observed for much larger action values. 

The splitting used in § was based on first splitting Cc 
and then l-L^ • Other engines can be obtained by different 
splitting of For example, consider the two-stroke 

engine obtained by splitting Cc + ^h' 

^2 stroke ^ ^-iU^c+U)^ ^ 

(14) 

Note that in the two-stroke engine the thermal coupling 
has to be | stronger compared to the continuous case 
in order to provide the same action. Using the SRT we 
obtain the complete equivalence relations of the three 
main engine types [77]: 


1^2 stroke 

stroke 


(15) 

|2 stroke ^ 
e(h) — 

/o4 stroke ^ 
^e(h) — 

/<4Cont 

^e(h)^ 

(16) 

j^2 stroke 

^4 stroke 

^cont 

(17) 


Another type of engine exits when the interaction with 
the work repository is carried out by two physically dis¬ 
tinct couplings. This happens naturally if 7 ^ 

E 2 — El so that two different driving lasers have to be 
used and the Hamiltonian is Hq + cos((F ^2 — Ei)t)Hwi + 
cos((F^4 — Es)t)Hyj 2 - In such cases, one can make the 
splitting shown in Fig. 5d. In this numerical example we 
use: H^i = e(t) |1) (2| + h.c. and 11^2 = e(t) |3) (4| + h.c. 
Since there are two different work strokes in addition to 
the thermal stroke this engine constitute a four-stroke 
engine which differece. 


C. Power and energy flow balance 


The average power and heat flow in the equivalence 
regime are independent of the cycle time: 


Pw = — = ( ifo 

Tr 


' eye 


2 ^“ 


Ps 


Je{h) 


%^ = (^o|£cw|p.>. 


(18) 

(19) 


Using the steady state definition © one obtains the 
steady state energy balance equation: 


Pw + Jc + — 0. 


( 20 ) 


Equation (20) does not necessarily hold if the system is 


not in steady state as energy may be temporarily stored 
in the baths or in the work repository. 
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small s large s 



Figure 7: Power as a function of action for various engine 
types in steady state. The 4 stroke variant (green) is described 
at the end of Sec. mE The action is increased by increasing 
the strokes duration (top illustration), (a) For large action 
with respect to h, the engines significantly differ in perfor¬ 
mance. In this example all engines have the same efficiency, 
but they extract different amounts of heat from the hot bath, 
(b) In the equivalence regime where the action is small, all 
engine types exhibit the same power and also the same heat 
flows. The condition s < h/2 that follows from the Strang de¬ 
composition agrees with the observed regime of equivalence. 
The time symmetric structure of the engines causes the devi¬ 
ation from equivalence to be quadratic in the action. 


Figure 7 shows the power in steady state as a function 
of the action. The action is increased by increasing the 
time duration of each stroke (see top illustration in Fig. 
7) . The field and the thermal coupling are e = jh = Ic = 
5 X 10“^. The coupling strengths to the bath and work 
repository are not changed. When the engine action is 
large compared to the engines behave very differently 
(Fig. 7a). On the other hand, in the equivalence regime, 
where s is small with respect to /i, the power of all engines 
types converges to the same value. In the equivalence 
regime the power rises quadratically with the action since 
the correction to the power is jrcyc ^ '^cyc' This power 
plateau in the equivalence regime is a manifestation of 
quantum interference effects (coherence in the density 
matrix), as will be further discussed in the next section. 

The behavior of different engines for large action with 
respect to h is very rich and strongly depends on the ratio 
between the field and the baths coupling strength. For 
example, if the field is amplified then for some parameter 
the four-stroke engine can produce more power than the 
continuous and two stroke engine. Some features of this 
diverse dynamics will be discussed elsewhere. 


Finally we comment that the same formalism and re¬ 
sults can be extended for the case the drive that is slightly 
detuned from the gap. 


D. Lasing condition via the equivalence to 
two-stroke engine 

Laser medium can thought of as continuous engine 
where the power output is light amplification. It is well 
known that lasing requires population inversion. Scovil 
et. al [55] were the first to show the relation between 
population inversion lasing condition, and the Carnot ef¬ 
ficiency. 

Using the equivalence principle, presented here, the 
most general form of the lasing condition can be obtained 
without any reference to light-matter interaction. 

Let us start by decomposing the continuous engine into 
an equivalent two-stroke engine. For simplicity, it is as¬ 
sumed that the hot and cold manifolds have some overlap 
so that in the absence of the driving field this bath leads 
the system to a unique steady state pQ. If the driving 
field is tiny with respect to the thermalization rates then 
the system will be very close to po in steady state. 

To see when po can be used for work extraction we 
need to discuss passive states. A passive state is a state 
that is diagonal in the energy basis, and with popula¬ 
tions that decreases monotonically with the energy |66) . 
The energy of a passive state cannot be decreased (or 
work cannot be extracted from the system) by applying 
some unitary transformation (the Hamiltonian after the 
transformation is the same as it was before the transfor¬ 
mation) [42j[66]. Thus, if po is passive, work cannot be 
extracted from the device regardless of the details of the 
driving field (as long as it is weak and the equivalence 
holds). 

A combination of thermal baths will lead to an energy 
diagonal po- Consequently, to enable work extraction, 
passivity must be broken by population inversion. There¬ 
fore, we obtain the standard population inversion condi¬ 
tion. Note that the derivation does not require Einstein 
rate equation and any information on the processes of 
emission and absorption of photons. 

Furthermore, it now becomes clear that if “coherent 
baths” are used uni so that Po is no longer diagonal in 
the energy basis (and therefore no longer passive) it is 
possible to extract work even without population inver¬ 
sion. 

In conclusion, using the equivalence principle it is pos¬ 
sible to import known results from work extraction in 
stroke schemes to continuous machines. 


IV. QUANTUM THERMODYNAMIC 
SIGNATURE 

Can thermodynamic measurement reveal quantum ef¬ 
fects in the engine? To answer this we first need to define 
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(b) 

|P)= 




f ' 



pop 


0 : 

h 

pop 


II 

II 

■ 

----- p. 



coh 


: 

0 

coh 

J 


'^pop: 

coh ^ 



Figure 8: (a left) Dephasing operations (slanted line, operator 
T>) commute with thermal baths so the dephased engine in (a) 
left, is equivalent to the one on the right. In the new engine 
the unitary evolution is replace by VIAV. If T> eliminates all 
coherences, the effect of T>IAT> on the populations can always 
be written as a doubly stochastic operator, (b) Any Hermitian 
Hamiltonian in Liouville space has the structure shown in (b). 
Thus, first order changes in populations critically depend on 
the existence of coherence. 


the corresponding classical engine. 

The term “classical” engine is rather ambiguous. There 
are different protocols of modifying the system so that is 
behaves classically. To make a fair comparison to the 
fully quantum engine we look for the minimal modifica¬ 
tion that satisfies the following conditions: 

1. The dynamics of the device should be fully de¬ 
scribed using population dynamics (no coherences, 
no entanglement). 

2. The modification should not alter the energy levels 
of the system, the couplings to the baths, and the 
coupling to the work repository. 

3. The modification should not introduce a new source 
of heat or work. 

To satisfy the first requirement we introduce a dephasing 
operator that eliminate the coherences!^ and leads to a 
stochastic description description of the engine. Clearly, 
a dephasing operators satisfy the second requirement. To 
satisfy the third requirement we require “pure dephasing”; 
a dephasing in the energy basis. The populations in the 
energy basis are invariant to this dephasing operation. 
Such a natural source of energy basis dephasing emerges 
if there is some scheduling noise m- That is, if there is 
some error in the switching time of the strokes. 

Let us define a “quantum-thermodynamic signature” 
as a signal that is impossible to produce by the corre¬ 
sponding classical engine as defined above. 

Our goal is to derive a threshold for power output that 
a stochastic engine cannot exceed but a coherent quan¬ 
tum engine can. 

Before analyzing the effect of decoherence it is instruc¬ 
tive to distinguish between two different work extraction 
mechanisms in stroke engines. 


A. Coherent and stochastic work extraction 
mechanisms 


Let us consider the work done in the work stroke of a 
two-stroke engine (as in Fig. 5c): 


IH = 





Writing the state as a sum of population and coherences 
1/5) = \ppop) + \hoh) we get: 


W = (^Ho 

+ (ho 


E 

n=l 

E 

n=l 




2n-l 


(2n- 1)! 


Pcoh 


{-i\Uwrw) 


2n 


(2n)! 


Ppop 


( 21 ) 


This result follows from the generic structure of Hamilto¬ 
nians in Liouville space. Any H that originates from an 
Hermitian Hamiltonian in Hilbert space (in contrast to 
Lindblad operators as a source) has the structure shown 
in Fig 8b (see Appendix II for Liouville space derivation 
of this property). That is, it connects only populations to 
coherences and vice versa, but it cannot connect popula¬ 
tions to populations directly [79]. In addition, since {Ho\ 
acts as a projection on population space one gets that 
odd powers of Hw can only operate on coherences and 
even powers can only operate on populations. Thus, the 
power can be extracted using two different mechanisms. 
A coherent mechanism that operates on coherences and 
a stochastic mechanism that operates on populations. 

The effect of the “stochastic” terms ^^=1 ^ — 

on the populations are equivalently described by a sin¬ 
gle doubly stochastic operator. If there are no coher¬ 
ences (next section) this leads to a simple interpretation 
in terms of full swap events that take place with some 
probability. 

Continuous engines, on the other hand, have only co¬ 
herent work extraction mechanism. This can be seen 
from the expression for their work output 


P 


cont 



1 


1 

2^- 

II 

2 ^“ 



( 22 ) 


where again we used the population projection property 
of (i7o|, and the structure of Pw (Fig- 8b). We conclude 
that in contrast to stroke engines, continuous engines 
have no stochastic work extraction mechanism. This dif¬ 
ference stems from the fact that the in steady state the 
state is stationary in continuous engines. Consequently, 
there are no higher order terms that can give rise to a 
population-population stochastic work extraction mech¬ 
anism. This is a fundamental difference between stroke 
engines and continuous engines. This effect is pronounce 
outside the equivalence regime where the stochastic terms 
become important (see Sec. [v| ). 
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B. Engines subjected to pure dephasing 


Consider the engine shown in Fig 8 a. The slanted 
lines on the baths indicate that there is an additional 
dephasing mechanism that takes place in parallel to the 
thermalization[ 8 Q] . Let us denote the evolution operator 
of the pure dephasing by V. In principle, to analyze the 
deviation from the coherent quantum engine, first the 
steady state has to be solved and then work and heat 
can be compared. Even for simple systems this is diffi¬ 
cult task. Hence, we shall take a different approach and 
derive a power upper bound for stochastic engines. It 
is important that the bound will contain only quantities 
that are unaffected by the level of coherence in the sys¬ 
tem. For example, average energy or dipole expectation 
values, do contain information on the coherence. We con¬ 
struct a bound in terms of the parameters of the system 
(e.g. the energy levels, coupling strengths, etc.), an is 
independent of the state of the system. In the pure de¬ 
phasing stage the energy does not change. Hence, the 
total energy change in the VUV stage is associated with 
work. 

Let Vcomp = \pop) {pop\ be a projection operator on 
the population space. This operator generates a complete 
dephasing that eliminates all coherences. In such case, 
the leading order in the work expression becomes 


W = 

^2 


'D p 

^comp^ 


comp 




(23) 


where we used {Ho\V = {Ho\ and Vcomp \p) = \Ppop)- 
Since Vcomp eliminates coherences, W does not contain 
a linear term in time. Next by using the following rela- 
tion: {Ho\B\p) < ^ (ifo |i?o) (p |/o) ||B||, ^ {Ho\h) = 
vt^) , we find that for s h the power of a stochastic 
engine satisfies: 


Pstoch < ytr{Hl)-tr{HoY^ld\cyc, (24) 
z = 1 two-stroke 

2 : = 1/2 four-stroke 


where is the gap of the interaction Hamiltonian (max¬ 
imal eigenvalue minus minimal eigenvalue of and 

d is the duty cycle - the fraction of time dedicated to 
work extraction (e.g. d = 1/3 in all the examples in 
this paper). We also used the fact that {ppop \ppop) is al¬ 
ways smaller than the purity {p\p) and therefore smaller 
than one. Note that, as we required, this bound is state- 
independent, and the right hand side of ( [^ contains no 
information on the coherences in the system. As shown 
earlier, in coherent quantum engines (in the equivalence 
regime) the work scales linearly with Tcyc (see ( p^ and 
(p!^) and therefore the power is constant as a function 
of '^cyc' When there are no coherences the power scales 
linearly with Tc^c- 


Numerical results of power as function of cycle time 
are shown in Fig. 9. The power is not plotted as func¬ 
tion of action as before, because at the same cycle time 
the coherent engine and dephased engine have different 
action. The action of the dephased engine is 


Sdeph — (||>^c|| + ||>C/i|| + 


\n. 


^^dephasing || )T 2 yc- (25) 


If the dephasing is significant the action is large and 
equivalence cannot be observed. That is, a fully stochas¬ 
tic engine in quantum system have large action and can¬ 
not satisfy s h. 

The chosen coherence time is lOOrd. When the cycle 
time is small with respect to the coherence time, equiva¬ 
lence is observed. Yet, the power is significantly smaller 
compared to the fully coherent case. For longer cycles 
the decoherence starts to take effect, and the expected 
linear power growth is observed. The stochastic power 
bounds for a two-stroke engine (dashed-blue), and for 
a four-stroke engine (dashed-red) define a power regime 
(shaded area) that is inaccessible to fully stochastic en¬ 
gines. Thus, any power measurement in this regime un¬ 
equivocally indicates the presence of quantum coherences 
in the engine. Note, that to measure power the work 
repository is measured and not the engine. Furthermore, 
the engine must operate for many cycles to reduce fluctu¬ 
ations in the accumulated work. To calculate the average 
power the accumulated work is divided by the total oper¬ 
ation time and compared to the stochastic power thresh¬ 
old @. 

Note that had we chosen complete dephasing then the 
power output of the continuous engine would have been 
zero as expected from ( 22 ). 

In summary, quantum thermodynamics signature in 
stroke engines can be observed in the weak action limit. 


V. THE OVER THERMALIZATION EFFECT IN 
COHERENT QUANTUM HEAT ENGINE 

In all the numerical examples studied so far, the uni¬ 
tary action and the thermal action were roughly com¬ 
parable for reason that will soon become clear. In this 
section we study some generic features that take place 
when the thermal action takes over. 

Let us now consider the case where the unitary con¬ 
tribution to the action is small with respect to 

h. All the time intervals are fixed but we can control 
the thermalization rate 7 (for simplicity, we assume it is 
the same value for both baths). Common sense suggests 
that increasing 7 should increase the power output. At 
some stage this increase will stop since the system will al¬ 
ready reach thermal equilibrium with the bath (or baths 
in two-stroke engines). Yet, Fig. 10 shows that there is 
a very distinctive peak where an optimal coupling takes 
place. That is, in some cases less thermalization leads 
to more power. We call this effect over thermalization. 
This effect is generic and not unique to the specific model 
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Figure 9: The output of the three types of engines (two-stroke 
blue, four stroke red, continuous black) with and without de¬ 
phasing (same as Fig. 7b). The dephasing time is lOOrdrive- 
Well above the dephasing time (~ 200Tdrive) the power grows 
linearly as expected from stochastic engines (three bottom 
solid lines). Below the dephasing time (~ ^Ordrive) equiva¬ 
lence is observed. Yet the power is significantly lower com¬ 
pared to the coherent engine (top three solid lines). The 
dashed lines show the stochastic upper bound on the power 
for two-stroke (dashed-blue) and four-stroke (dashed-red) en¬ 
gines. Any power measurement in the shaded area of each 
engine indicates the presence of quantum interference in the 
engine. This plot also demonstrates that for weak couplings 
(low action) coherent engines produce much more power com¬ 
pared to stochastic dephased engines. 


used in the numerical simulations. The parameter used 
for the plot are e = 7 c = 7 /i = 2 x 10 “^ and the number 
of drive cycle pen engine cycle is m = 600. 

The peak is a consequence of the interplay between the 
two different work extraction mechanisms (see Sec. 4.1). 
For low 7 the coherences in the system are significant 
and the leading term in the power is (F7o | \ pcoh) d 

(where d is the duty cycle). In principle, all Lindblad 
thermalization processes are associated with some level 
of decoherence. This decoherence generates an expo¬ 
nential decay of \pcoh) that explains the decay on the 
right hand side of the peak. At a certain stage the lin¬ 
ear term becomes so small that the stochastic second or¬ 
der term — | (^Hq \ ppop) dP^cyc dominates the power. 
\ppop) eventually saturates for large 7 and therefore the 
stochastic second order term leads to a power saturation. 
Interestingly, in the example shown in Fig. 0 we ob¬ 
serve that the peak is obtained when 7 and e are roughly 
equal. Of course, what really matters is the thermal ac¬ 
tion with respect to unitary action and not just the values 
of the parameter 7 and e. 

If thermalization occurs faster, the thermal stroke can 
be shortened and this increases the power. However, 
this effect is small with respect to the exponential de¬ 
cay of the coherences. We conclude that even without 
additional dephasing as in the previous section, exces¬ 
sive thermal coupling turns the engine into a stochastic 
machine. For small unitary action this effect severely de¬ 



Figure 10: The over thermalization effect is the decrease of 
power when the thermalization rate is increased. Over ther¬ 
malization degrades the coherent work extraction mechanism 
without effecting the stochastic work extraction mechanism. 
When the coherent mechanism gets weak enough, the power 
is dominated by the stochastic power extraction mechanisms 
and power saturation is observed (dashed lines). The contin¬ 
uous engine has no stochastic work extraction mechanism and 
therefore it decays to zero without reaching saturation. 

grades the power output. The arguments presented here 
are valid for any small action coherent quantum engine. 


VI. CONCLUDING REMARKS 

We identified coherent and stochastic work extraction 
mechanisms in quantum heat engines. While stroke en¬ 
gines have both, continuous engines only have the coher¬ 
ent mechanism. We introduced the “norm action” of the 
engine using Liouville space and showed that when this 
action is small compared to h all three engine types are 
equivalent. This equivalence emerges because for small 
action only the coherent mechanism is important. De¬ 
spite the equivalence, before the engine cycle is completed 
the state of the different engine type differ by O(s^). This 
holds true also for work and heat. Remarkably, at the 
end of each engine cycle a much more accurate 0 ( 5 ^) 
equivalence emerges. Furthermore, the equivalence holds 
also for transient dynamics, even when the initial state is 
very far from the steady state of the engine. It was shown 
that for small action the coherent work extraction is con¬ 
siderably stronger than the stochastic work extraction 
mechanism. This enabled us to derive a power bound for 
stochastic engines that constitutes a quantum thermody¬ 
namics signature. Any power measurement that exceeds 
this bound indicated the presence of quantum coherence 
and the operation of the coherent work extraction mech¬ 
anism. 

The present derivation makes no assumption on the 
direction of heat flows and the sign of work. Thus our 
result are equally applicable to refrigerators and heaters. 

It is interesting to try and apply these concepts of 
equivalence and quantum thermodynamic signatures to 
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more general scenarios: non Markovian baths, engines 
with non symmetric unit cell, and engines with correla¬ 
tion between different particles (entanglement and quan¬ 
tum discord). We conjecture that in multiple particle 
engines entanglement will play a similar role to that of 
coherence in single particle engines. 

Work support by the Israeli science foundation. Part 
of this work was supported by the COST Action MP1209 
’Thermodynamics in the quantum regime’. 


Appendix I - Strang decomposition validity 

Let JC be an operator generated by two non commuting 
operators A and B: 


The splitted operator is 

_ ^^Adt^Bdt^^Adt 


Our goal is to quantify the difference between 1C and 
/Cs, ||/Cs —/C|| where H-H stands for the spectral norm. 
In principle, other sub-multiplicative matrix norm can 
be used (such as the Hilbert Schmidt norm). However, 
the spectral norm captures more accurately aspects of 
quantum dynamics EHHZI]- 1C can be expanded as: 

K-s on the other hand is; 




E 


B^dt^ {Al2)^dt^ 

k\ T\ m\ 


oo n n—l 


EEE 


{Al2f & 

k\ l\{n-l-k)\ 


de. (27) 


Due to the symmetric splitting the terms up to n = 2 (in¬ 
cluding n=2) are identical for both operators. Therefore 
the difference can be written as 


||/C,-/C|| 


OO n n—l 

EEE 


n=3 1=0k=0 


{A/2f B^ {Al2r-^-^ 

k\ l\ {n-l-k)\ 


(M H“ B^^dt' 

^ n\ 

n=3 


next we apply the triangle inequality and the sub- 
multiplicativity property and get: 


||/C,-/C|| < 

^^^ \\A/2\\>^\\B\t\\A/2\r-^-’^ 
fc! l\ {n-l-kV. 

n=3 1=0 k=0 ' ^ 

, ^(P|| + ||B||)"df" 

+ T-m- • (29 


Using the binomial formula two times one finds 

VvC 11-^/211 "ll^ll' IW2|r~'~" ,.n ^ 

Z^Z^Z^ fc! n in-l-k)\ 

n=3 1=0 k=0 ' ' 

g {\\A\\ + \\B\\rdt- ^ 

n=3 

and therefore 

Wfc^-icw < 

n=3 

= 2RMA\\ + \\B\\)dt]. (31) 

The right hand side is the Taylor reminder of a power 
series of an exponential with as argument s = {\\A\\ + 
||S||)(it. The Taylor reminder formula for the exponent 

function is Rk{x) = (feVi)! where 0 < ^ < I (for now 
we assume s < 1). Setting k = 2 and ^ = I (worst case), 
we finally obtain 

\\ICs-m < ^[{\\A\\ + \\B\\)dtf<s^ (32) 

s = (||A|| + ||B||)df. (33) 

To get an estimation where the leading non neglected 
term of /C, {A + S)^dt^/2 is larger then the reminder we 
require that 


{\\A^B\\fdty2 > s\ (34) 

Using the triangle inequality we get the estimated condi¬ 
tion for the Strang decomposition: 

5 < 1/2. (35) 

This condition explains why it was legitimate to limit the 
range of s to I in the reminder formula. 


Appendix II - Liouville space formulation of 
quantum dynamics 

Quantum dynamics is traditionally described in 
Hilbert space. However, it is convenient, in particular for 
open quantum systems, to introduce an extended space 
where density operators are vectors and time evolution 
is generated by a Schrodinger-like equation. This space 
is usually referred to as Liouville space [62]. We denote 
the "density vector" by |p) G . It is obtained by 

reshaping the density matrix p into a larger single vector 
with index a G {1,2, The one-to-one mapping of 

the two matrix indices into a single vector index {i, j} ^ 
a is arbitrary, but has to be used consistently. The vec¬ 
tor Ip) is not normalized to unity in general. Its norm is 
equal to the purity, V = tr(p^) = (p |p) where (p| = |p)^ 
as usual. The equation of motion of the density vector in 
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Liouville space follows from dtpa = ppd{dtpa)ldp^. 
Using this equation, one can verify that the dynamics of 
the density vector |r) is governed by a Schrodinger-like 
equation in the new space, 


idt \p) ='H\p), (36) 

where the super-Hamiltonian G ^ is given by, 


n. 


aj3 


.d{dtpa) 

' dpp ■ 


(37) 


A particularly useful index mapping is described in m 
and in [65]. For this form 1-L has a very simple form in 
terms of original Hilbert space Hamiltonian and Lindblad 
operators. 

1-L = li^ + £ is non-Hermitian for open quantum sys¬ 
tems. li^ originates from the Hilbert space Hamiltonian 
i7, and C from the Lindblad terms. is always Hermi- 
tian. The skew-Hermitian part (£ — £’t')/2 is responsible 
for purity changes. Yet, in Liouville space, the Lindblad 
operators A]^ 0 may also generate an Hermitian term 
(£ + Though Hermitian in Liouville space this 

term cannot be associated with a Hamiltonian in Hilbert 
space. 

For time-independent H the evolution operator in Li¬ 
ouville space is: 


\p{t))=KW))=e-^^^^-^')\p{t')). (38) 


If £ = 0, /C is unitaryThe fact that the evolution opera¬ 
tor can be written as an exponent of a matrix, without 
any commutators as in Hilbert space, is a very significant 
advantage (see for example [68]). It is important to note 
that not all vectors in Liouville space can be populated 
exclusively. This is due to the fact that only positive 
p with unit trace are legitimate density matrices. The 
states that can be populated exclusively describe steady 
states, while others correspond to transient changes. We 
remind that in this paper we will use calligraphic letters 
to describe operators in Liouville space and ordinary let¬ 
ters for operators in Hilbert space. For states, however, 
\A) will denote a vector in Liouville space formed from 
^NxN by “vec-ing” A into a column in the same proce¬ 
dure p is converted into \p). 


The proof is as follows: 


njH 

^ij,mn 


Himdjn 




(40) 


Therefore, using we get: 

u^\h) = Tu%h^ = Th^ 


mn 




(41) 


This property is highly useful. We stress that (39) is a 


property of Hermitian operators in Hilbert space where 
both H and H are well defined. A general Hermitian 
operator in Liouville space may not have a corresponding 
H in Hilbert space. 


Another property that immediately follows from (40) 


IS 


^Ikk = 0. (42) 

This corresponds to a well known property of unitary 
operation. If the system starts from a diagonal density 
matrix, then for short times the evolution generated by 
= / — iH^dt + 0{dt^) does not change the 
population in the leading order. 


Expectation values and their time evolution in 
Liouville space 


The expectation value of an operator in Hilbert space is 
{A) = tr{pA) . Since p is Hermitian the expectation value 
is equal to the inner product of A and p and therefore: 


(^) =tr{pA) = {p\A). 

the dynamics of {A) under Lindblad evolution operator: 


^^{A) = -i{A\n\p) + ^P 


u 

dt 


(43) 


Note that in Liouville spaee there is no eommutator term 
since H operates on \p) just from the left. If the total 
Hamiltonian is Hermitian and time independent the con¬ 
servation of energy follows immediately from applying 
(|4^ and (ra for A = H. 


Useful relations in Liouville space. 


Appendix III - The symmetric rearrangement 
theorem (SRT) 


In Liouville space, the standard inner product of two 
operators in Hilbert space tiA^B reads 

tvA'<B= {A\B). 

In particular the purity V = (r |r) is just the square of 
the distance from the origin in Liouville space. 

A useful relation for : 

(39) 


The goal of this Appendix is to explain why the equiv¬ 
alence of evolution operators leads to equivalence of work 
and equivalence of heat. In addition it is shown why it is 
valid also for transients. For the equivalence of evolution 
operator we have required that the super Hamiltonian is 
symmetric and that the action is small: 

n{t) = n-t), (44) 

i-+r/2 

s = I \\1-L\\dt<^h. (45) 

i-r/2 


H^\H) = {H\H^ = Q. 
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Let the initial state at time t = — t/2 be 


II 

(46) 

this state leads to a final state at r/2 


\h) = lp('r/2) ■ 

(47) 


Our goal is to evaluate a symmetric expectation value 
difference of the form: 


dAfot — [(^(^ 2 )) - (^(^ 1 ))] + [(^(-^ 1 )) - (^(-^ 2 ))] 

= [{A\p{t,)) - {A\p{t,))] 

+ [(A|p(-ti))-(A|p(-t2))], (48) 

t2,ti > 0 


that is, the change in the expectation value of A in the 
segment [^ 1 ,^ 2 ] and its symmetric counterpart in negative 
time (e.g. the green areas in Fig. 11a). When A is equal 
to Hq this difference will translate into work or heat. We 
start with the expansion: 


[(^(^ 2 )) - (^(^ 1 ))] — {A “ d\p{ti)) — 

A —i'H{ti)6t — -'H{ti)‘^St‘^\p{ti) \ + 0{s^). 


(49) 



Figure 11: The Hamiltonians in (a) and (b) are related by 
symmetric rearrangement of the time segments. Up to a small 
correction O(s^) the change in expectation values of an ob¬ 
servable A that takes place during the green segments is the 
same in both cases. This effect explains why work and heat 
are the same in various types of engines when s is small com¬ 
pared to h (equivalence regime). 


For the negative side we get: 

liM-ti)) - {A{-t 2 ))] = {A\I- \r{-ti)) = 

A ^ + O(s^). 

(50) 

Next we use the fact that: 

\p{ti)) = \K0))-i[ 'H{t)dt\p{0))+0{s^), 

Jo 

(51) 

\K-ti)) = \p{o)) + i[ niOdtipio))+ o{s^). 

Jo 


(52) 


When adding the two segments the second order cancels 
out and we get: 

SAtot = -2i {A \H{ti)\p{Q)) 5t + 0{s^). (53) 


We use 1-L to denote a Liouville space operator 
which may be any unitary operation or Markovian Lind- 
blad operation. Assume that in addition to the symmet¬ 
ric bins of interest (green bins) the reminder of the time 
is also divided into bins in a symmetric way so that there 
is still a reffection symmetry also in the bin partitioning. 
Now, we permute the bins in the positive side as desired, 
and then make the opposite order in the negative side so 
that the reffection symmetry is kept. An example of such 
an operation is shown in Fig. 11b. Due to the Strang de¬ 
composition we know that the total evolution operator 
will stay the same under this rearrangement up to third 
order: 

+ O(s^), (54) 

where Tsym [x] stands for evaluation of x after a symmet¬ 
ric reordering. 


The symmetric rearrangement theorem (SRT) 


Note that the result using |p(0)) which is not given ex¬ 
plicitly. To correctly relate it to |p(— t/ 2)) we have to 
use symmetric rearrangement properties of the evolution 
operator. 


Symmetric rearrangement 

In Fig. 11a there is an illustration of some time de¬ 
pendent Hamiltonian with reffection symmetry H(t) = 


From ( [ 5 ^ we get that if the initial state is the same for 
a system described by /C, and for a system described by 
Tsym[JQ^ the final state at t = r/2 is the same for both 
systems up to a third order correction: 


Kl))=Tsy^[\Kl))] + Ois% (55) 


Using (51 ),(|52[) we see that: 


I-(o)j, 


(56) 
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and because of (55) it also holds that: 

Tsym[\p{^))] = |p(0))+O(5^) = l^(-^ 


+ 0 (s^), 

(57) 


using this in (53) we get that: 

6Atot = -2i {A\n{ti) 




St + Ois-^). (58) 


Expression (58) no longer depends on the position of the 
time segment, but only on its duration and on the value 
of H. Thus, the SRT states that the expression above 
also holds for any symmetric rearrangement 


dAiQi — 'T'sym\dAiQ-^ + 0(s ). 


(59) 


If we replace A by Hq and 'H(ti) by Cc^Ch or 1-Lw 
we immediately get the invariance of heat and work to 
symmetric rearrangement (up to 5^). If |p(—^)) is the 
same for all engines then |p(^)) is also the same for all 
engines types up to 0{s^). Consequently for all stroke 
engines the expression for work and heat are: 

W = -2i{A\ J +0(s3), 

tet^ 


(60) 

Qc{h) = -2i(A| J Ap( 2 )) 

(61) 


Using the identity |p(i)H |p(-i)) = |p(0) + |p(-0) + 
O(s^) that follows from ( [^ , the integration over time 
of the energy flows = (i?o \ ^'Hw\Kt)) and jc(h) = 
(Hq \ £c(h) I p{t)) for continuous engines, yields expres- 
sions ([6Q| and (|M] ) once more. This implies that the 
SRT ( |60p and (|6ip holds even if the different operations 
>Cc, jCh and Hw overlap with each other. 


We emphasize that all the above relations hold for 
any initial state and not only in steady state where 
|p(|)) = |p(—D). The physical implication is that in 
the equivalence regime different engines are thermody¬ 
namically indistinguishable when monitored at the end 
of each cycle, even when the system is not in its steady 
state. 
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